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Understanding and controlling the shape of thin, soft objects has been the focus of significant 
research efforts among physicists, biologists, and engineers in the last decade. These studies aim to 
utilize advanced materials in novel, adaptive ways such as fabricating smart actuators or mimicking 
living tissues. Here, we present the controlled growth-like morphing of 2D sheets into 3D shapes 
by preparing geometric composite structures that deform by residual swelling. The morphing of 
these geometric composites is dictated by both swelling and geometry, with diffusion controlling the 
swelling-induced actuation, and geometric confinement dictating the structure’s deformed shape. 
Building on a simple mechanical analog, we present an analytical model that quantitatively de¬ 
scribes how the Gaussian and mean curvatures of a thin disk are affected by the interplay among 
geometry, mechanics, and swelling. This model is in excellent agreement with our experiments and 
numerics. We show that the dynamics of residual swelling is dictated by a competition between 
two characteristic diffusive length scales governed by geometry. Our results provide the first 2D 
analog of Timoshenko’s classical formula for the thermal bending of bimetallic beams - our gener¬ 
alization explains how the Gaussian curvature of a 2D geometric composite is affected by geometry 
and elasticity. The understanding conferred by these results suggests that the controlled shaping of 
geometric composites may provide a simple complement to traditional manufacturing techniques. 

PACS numbers: 


The continuous shape change during the growth and 
decay of biological structures is a constant presence 
within the natural world. These structures morph to 
accommodate an influx of new material, either growing 
due to an increase of mass to the system from an exter¬ 
nal source, or swelling from the absorption of excess fluid, 
such as water, caused by a change in humidity. Often, 
these morphological changes from swelling or growth re¬ 
sult in geometries that enhance the biological structure’s 
functionality m, for instance the Venus flytrap’s leaves 
snap closed after osmotically swelling, and this structural 
reconfiguration is essential for its nutrition [3]. Some 
of the most dramatic growth-induced deformations oc¬ 
cur with slender structures, such as growing leaves [4], 
wrinkling skin 0, and the writhing of tendril-bearing 
climbers [6]. Thin structures like these significantly de¬ 
form to adopt nontrivial three dimensional shapes be¬ 
cause they must bend to release their stretching en¬ 
ergy 0 . The coupling between growth and large de¬ 
formations presents an interesting opportunity for the 
morphing of synthetic structures, whereby if specific re¬ 
gions within a thin material can be prescribed to stretch 
the overall structure will adopt a new shape. Here, we 
demonstrate the controlled morphing of a thin structure 
by locally altering its intrinsic geometry; introducing a 
new class of geometric composite structures. While com¬ 
posite materials combine constituents to enhance physi¬ 
cal or chemical properties, geometric composites combine 
different intrinsic geometries to produce shapes that dif¬ 
fer from the individual components. Shape change in 
these geometric composites is triggered by the growth¬ 
like swelling of one region due to residual fluid in the 
surrounding material. By tuning the underlying materi¬ 


als and geometry, we present a novel and straightforward 
means to directly morph elastic materials into various 3D 
shapes. This growth-actuated manufacturing provides a 
natural complement to traditional reductive [8] and ad¬ 
ditive manufacturing techniques isi mg. 

The continuous shape change during the growth and 
decay of biological structures is a constant presence 
within the natural world. These structures morph to 
accommodate an influx of new material, either growing 
due to an increase of mass to the system from an exter¬ 
nal source, or swelling from the absorption of excess fluid, 
such as water, caused by a change in humidity. Often, 
these morphological changes from swelling or growth re¬ 
sult in geometries that enhance the biological structure’s 
functionality m, for instance the Venus flytrap’s leaves 
snap closed after osmotically swelling, and this structural 
reconfiguration is essential for its nutrition |3]. Some 
of the most dramatic growth-induced deformations oc¬ 
cur with slender structures, such as growing leaves [4], 
wrinkling skin [5], and the writhing of tendril-bearing 
climbers [6]. Thin structures like these significantly de¬ 
form to adopt nontrivial three dimensional shapes be¬ 
cause they must bend to release their stretching en¬ 
ergy 0 - The coupling between growth and large defor¬ 
mations presents an interesting opportunity for the mor¬ 
phing of synthetic structures, whereby if specific regions 
within a thin material can be prescribed to stretch the 
overall structure will adopt a new shape. Here, we show 
the controlled morphing of a thin structure by locally al¬ 
tering its intrinsic geometry; preparing structures that 
are effectively geometric composites. While composite 
materials combine constituents to enhance physical or 
chemical properties, geometric composites combine dif- 
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ferent intrinsic geometries to produce shapes that differ 
from the individual components [71 fTTj . Shape change in 
these geometric composites is triggered by the growth- 
like swelling of one region due to residual fluid in the 
surrounding material. By tuning the underlying materi¬ 
als and geometry, we present a novel and straightforward 
means to directly morph elastic materials into various 3D 
shapes. This growth-actuated manufacturing provides a 
natural complement to traditional reductive [8] and ad¬ 
ditive manufacturing techniques 0110 ]. 

Understanding the interplay among geometry, mechan¬ 
ics, and swelling is essential for the controlled morph¬ 
ing of thin structures. Geometrically, Gauss’s Theorema 
Egregium states that the Gaussian curvature of a sheet is 
dictated by the local distances between points, and this 
intrinsic geometry is conserved by isometries. Therefore, 
if a stimulus changes these local distances, thereby chang¬ 
ing the sheet’s Gaussian curvature, the disk will most 
likely bend into a three dimensional shape. The under¬ 
lying geometry of a sheet may be programmed mechani¬ 
cally through confinement, or it can be altered chemically 
by locally swelling regions of the material. For swelling 
to induce an out-of-plane deformation it must stretch 
the sheet in a non-homogenous manner, either via selec¬ 
tive surface wetting m or anisotropic material proper¬ 
ties m- To describe the shape of the resulting deformed 
structures, the Foppl-von Karman equations were gener¬ 
alized to account for growth [13], and employed to study 
the morphogenesis of a blooming lily mi, leaves m and 
the buckling of swelling gels m ■ Building on these ideas 
of elasticity and growth, the theory of non-Euclidean 
plates mm was developed without any a priori as¬ 
sumptions regarding the structure’s displacements. This 
approach provided a powerful way to design and study 
the morphing of thin structures such as the buckling of 
sheets [18j Iff and ribbons [20] and the three dimensional 
transformations of hydrogels into helices m- Dias et al. 
presented an elegant study of the inverse buckling prob¬ 
lem, where the growth patterns corresponding to some 
desired axisymmetric shapes were found explicitly [22] . 
In both models, an analytical prediction of a morphing 
structure’s shape is often elusive, even for simple geome¬ 
tries and symmetric shapes. These concepts of differ¬ 
ential growth have been investigated in the ID case of 
morphoelastic rods, where instabilities such as the cir¬ 
cumferential buckling of a growing cylinder occur [23] . 
This work will build upon the contributions from differ¬ 
ential geometry and elasticity to analytically predict the 
shape of 2D geometric composites that change shape in 
response to residual swelling. 

In Fig. [l] we show a geometric composite that morphs 
from a flat plate into a negatively curved surface from the 
residual swelling of the annulus (green) by free polymer 
chains present in the disk (pink). To predict this struc¬ 
ture’s shape change, we first constructed a mechanical 
analog based on the geometry of the geometric compos¬ 
ite shown in figure [l] In our mechanical analog, the sheet 
is made of a disk and an annulus that share the same ma- 
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FIG. 1: (a) An initially flat, elastic plate morphs into a neg¬ 
atively curved shell, or saddle, in response to swelling of the 
annulus from residual free chains present in the inner disk, 
(b) A schematic representation of 2D plate morphing into a 
3D shell. 


terial properties, but are geometrically incompatible: the 
inner radius of the annulus is slightly bigger or smaller 
than the radius of the disk so that the annulus has to be 
mechanically strained for the composite to be bonded to¬ 
gether. The resulting geometrical frustration caused the 
disk to bend into positively or negatively curved surfaces, 
that is, into surfaces with positive or negative Gaussian 
curvature, respectively. This mechanical analog enforces 
the importance of geometry, and poses a common prob¬ 
lem in elasticity usually cited as a prototypical example 
for frustrated structures, but has never been fully inves¬ 
tigated [24] . In this work, we present an analytical model 
based on a membrane approximation that quantitatively 
captures the Gaussian curvature of the disks, and under¬ 
scores the strong connection between the Gaussian cur¬ 
vature and the metric. The agreement among numerics, 
experiments, and analytics is excellent. Consequently, we 
apply our model to explain residual swelling, and study 
how the Gaussian curvature of the disk varies with the ge¬ 
ometry of the composite. We show how the coupling be¬ 
tween swelling and geometry dictates the main feature of 
structural morphing, and we provide the first 2D exten¬ 
sion of Timoshenko’s formula for the thermal bending of 
bimetallic beams [25]. Finally, we describe the dynamics 
of residual swelling, which exhibits a transition between 
two characteristic diffusive length scales depending on 
the geometry of the composite. 
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Intrinsic &; Extrinsic Geometry 

In thin structures, the effect of mechanical or chemical 
stimuli can be represented mathematically by a target 
metric a - a tensor that describes what the local dis¬ 
tances between points should be in response to the ap¬ 
plied stimuli [IT. The actual distances between points 
across a surface are represented by the realized metric 
a and uniquely determine its Gaussian curvature, which 
describes the product of the surface’s two principal cur¬ 
vatures, k i and Hv 2 • If the Gaussian curvature imposed 
by the stimuli is not zero, the target metric cannot be 
realized in a flat disk, and the disk will bend to mini¬ 
mize the in-plane strains. Since the ratio of the bending 
to stretching energies scales like Ub/Uj ~ h 2 , the disk 
will bend to minimize the in-plane strains if the thick¬ 
ness h is small 1 . Therefore, for thin disks, a reasonable 
approximation is to determine the structure’s shape by 
minimizing the stretching energy, and allowing the disk 
to bend in whatever way necessary to accommodate the 
target metric. We write the bending energy density as 
Eh 3 (AH 2 — K), where H = (n\ + /c 2 )/2 is the mean cur¬ 
vature of the mid-surface of the sheet, an extrinsic ge¬ 
ometry property that depends on the embedding space, 
K — k\K 2 is the Gaussian curvature, an intrinsic geomet¬ 
ric property of the surface, and E is the Young modulus 
of the sheet. Then, considering an incompressible ma¬ 
terial (v = 0.5), the stretching energy may be written 
as [17] 

U s ~h E[tr 2 (a — a) + tr (a — a) 2 } y/\&\ dA . (1) 


I. GEOMETRIC COMPOSITES 
A. Membrane Approximation 

Before we consider the dynamic problem of a swelling- 
induced changing metric, we construct as a proof of con¬ 
cept a mechanical analog of the geometric composite and 
consider the static case of a mechanically prescribed hy¬ 
perbolic or elliptic metric over a disk. In the simplest 
case, when homothetic transformation radially stretches 
a circular disk, the radial and azimuthal distances all 
expand uniformly, and the disk stays flat since the sur¬ 
face’s metric remains Euclidean (Fig.[2^). If we consider, 
instead, a circular disk of radius R and an annulus of in¬ 
ner radius Ri and outer radius R e (Fig.J^), we can then 
impose a homothetic transformation independently on ei¬ 
ther the disk or the annulus. By defining a = Ri/R as 
the mismatch between the disk and the annulus, it is im¬ 
mediately apparent that if a ^ 1 the two structures are 


incompatible and the annulus must be stretched (a < 1) 
or compressed (a > 1) to fit the circular disk (Fig. [2^) . 
When the inner radius of the annulus is bonded to the 
edge of the circular disk, the resulting disk is a geomet¬ 
ric composite that may be roughly modeled as a body 
having the following target metric in polar coordinates 


a = /2 ( r ) (o r 2 ) ’ 


1, r < R 
a 2 , r > R 


( 2 ) 


This metric is flat within each of the two domains, but the 
metric of the geometric composite is not flat, i.e. there 
does not exist any parabola that fits aoo(r ), the azimuthal 
covariant metric coefficient. As in[26l we approximate all 
the strains to zero but ago — dgg; therefore, if the disk and 
the annulus are made of the same material, the stretching 
energy from equation 0 reads (see Appendix) 


>/ 




— a 2 r 2 ) 2 


dr . 


( 3 ) 

Physical intuition tells that when the annulus is 
stretched (compressed), the disk will bend into a dome¬ 
like (saddle-like) shape. This statement may be math¬ 
ematically represented as sgnA = sgn(l — a). To 
describe the resulting shape, we use Gaussian normal 
coordinates (p, 6) to express the realized metric [27] . 
where p(r) = J Q r yja rr (r')dr' measures the arc length 
along radial geodesics while 0 is the azimuthal angle. 
In these coordinates, the first fundamental form is writ¬ 
ten as ds 2 = dp 2 + CLQo(p)d6 2 , and by the Gauss theo¬ 
rem, the Gaussian curvature is —d pp sjagg / y/agg, where 
d pp is the second order partial derivative with respect 
to p. We minimize the stretching energy by looking 
for metrics with constant Gaussian curvature, that is 
aoo(p) = (sin (y/K p) / y/K) 2 . 2 As long as \K\ < 1 /R 2 e 3 , 
we can Taylor expand age (p) to linearize the metric in K 
as 


aee{p) = p 2 ~ y / + 0(^ 5 ) • (4) 

Note that the first order term corresponds to a flat met¬ 
ric whereas the second one dictates the kind of non- 
Euclidean geometry that the disk will develop depend¬ 
ing on the sign of K. The energy is quadratic in K 
and therefore can be minimized analytically; notice that, 
if the annulus is neither deformed (a = 1) nor present 
(Ri = R e ), the energy is a simple parabola in K with 
the minimum at K = 0 since the disk is not constrained, 
and does not need to bend. Similarly, when the disk is 
not considered (R = 0), the annulus does not need to 


1 The elastic energy is U = hU s +h 3 Ui) where U s has the dimension 

of an energy over a length, while lAb has the dimension of an 

energy over the cube of a length. 


2 Notice that when the Gaussian curvature is nega¬ 
tive, i.e. K < 0, the metric may be rewritten as 

aoo(p) = (sinh p) / y/^K) 2 . 

3 This upper bound means that each principal direction cannot 
have a curvature that exceeds 1/R e • 
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bend either, and stays flat (K = 0) with a radial stretch 
equal to a. Once the stretching energy is minimized, 
we observe that the bending energy density is equal to 
3 H 2 for dome-like shapes (K = H 2 ) and to 4 H 2 — K 
for saddle-like shapes. In the latter case, since K < 0, 
the disk tries to morph into a minimal surface (H = 0). 
Important aspects, albeit beyond the scope of this work, 
are the study of the discontinuities in the metric and the 
effect of a finite thickness. 

When the target metric is elliptic, the resulting shape 
is unique m • On the other hand, when the target metric 
is hyperbolic, the embedding is not unique, and shapes 
that are more complex then a saddle may develop when 
the thickness is very small [28] . In our case, both exper¬ 
imental and numerical evidence indicate that the thick¬ 
ness to radius ratio (~ 0.16) is sufficiently high to avoid 
the development of complex shapes other than the sad¬ 
dle, yet small enough for the structure to be considered 
thin. 


II. MECHANICAL ANALOG 
A. Experiments and Numerics 

To test our model, we prepared geometrically frus¬ 
trated structures to realize dome-like and saddle-like 
disks, and measured their Gaussian curvatures. Circu¬ 
lar molds were laser cut out of acrylic sheets, and used 
to cast samples with polyvinylsiloxane (PVS - Zhermack 
Elite Double 32). In these model experiments, we use 
an elastomer with a Young modulus E = 0.96 MPa and 
a Poisson ratio v = 0.5. The molds had a thickness 
h = 1.6 mm with the radii of the disks varying between 5 
and 12 mm. We designed the geometry so that the outer 
radius of the stretched disks was equal to 10 mm. Bond¬ 
ing between the stretched annulus and the disk was ac¬ 
complished by a small amount of uncrosslinked PVS. To 
measure their Gaussian curvature, we projected a laser 
sheet normal to the disk, and captured images of the re¬ 
flected light with a Edmund-Optics GigE camera with a 
Nikkor lens (35 mm f: 1-1.4) at 24 equally spaced points 
along the disk’s diameter. Image analysis was performed 
using Matlab to reconstruct the deformed shape. The 
annuli of polyvinylsiloxane elastomer were stretched ho- 
mothetically and bonded to the inner disk. Upon re¬ 
lease from the molds, the disks spontaneously morphed 
into domes or saddles. The annuli may be thought of 
as springs that want to release the energy by recover¬ 
ing their original shapes: for example, figure [2] (c) shows 
how the annulus must be compressed to get a saddle. 
We measured the shape of the deformed disks to deter¬ 
mine the two principal radii of curvature, and plotted 
the Gaussian curvature in figure |2j These measurements 
confirmed the assumption used in our analysis that K is 
approximately constant across the disk. We also carried 
out numerical simulations to solve the problem within 
the context of finite incompressible tridimensional elas¬ 


ticity with large distortions using a Neo-Hookean mate¬ 
rial model [29 implemented in the commercial software 
COMSOL Multiphysics. The stimulus in the numeri¬ 
cal model is represented by a unimodular distortion field 
F 0 = /(r)(I— e 3 ®e 3 )+f(r)~ 2 e 3 ®e 3 , where e 3 is the unit 
vector field orthogonal to the undeformed mid-surface of 
the disk. In this way, the tridimensional model is consis¬ 
tent with the 2D model since F^F 0 • e 7 0 = a 7M , 

where e 7 is the y-th vector of the covariant basis span¬ 
ning the undeformed mid-surface of the disk. The nu¬ 
merical results are also plotted in figure [2] 


B. Analysis 

The main plot in figure [2] refers to the case when 
R/R e a = 0.7. The blue curves show the azimuthal com¬ 
ponent of the target metric of the mid-surface agg as a 
function of r. As long as r < R, age is equal to p 2 since 
the inner disk has not been stretched; on the contrary, 
within the annulus (r > /£), the azimuthal component 
is equal to a 2 p 2 with a < 1 (a > 1) if the annulus has 
been stretched (compressed). Notice that the parabola 
in the inner disk is represented also for r > R (dashed 
blue curve) to show how the metric should look like if 
the disk were flat and not stretched. By the Taylor ap¬ 
proximation of the metric, it is evident that if K > 0 
(K < 0), agg stays below (above) the parabola. Solid 
black lines are the analytical solutions of Q for ellip¬ 
tic and hyperbolic target metrics. By looking at equa¬ 
tion ©, we notice that for the energy to be minimized, 
the analytical solution should be as close as possible (in 
a L 2 sense 4 ) to the target metric with weight functions 
that are r -3 and a~ 2 r ~ 3 in the disk and in the annulus, 
respectively. This explains why the realized metric stays 
in between the target metrics, and closer to the target 
metric in the disk. Circles and triangles in the main plot 
show the experimental result (we measure K and com¬ 
pute agg from equation Q). Dashed red curves represent 
the numerical solutions and show that the assumption of 
homogenous Gaussian curvature accurately describes the 
metric in this case apart from slight deviations near the 
edges. The experimental and numerical results are in 
excellent agreement with our closed form analytical solu¬ 
tion. The numerics also shows that H ~ \/K for domes 
and H ~ 0 for saddles, as predicted analytically. It is 
interesting to note that the capability of estimating an 
extrinsic geometric quantity, i.e. the mean curvature i7, 
from a model built around intrinsic geometries allows us 
to determine the displacement field, i.e. the embedding 
of the structure, up to rigid motions. This indetermi¬ 
nation does not affect the predictability of axisymmetric 
shapes like domes, but it does affect the one of saddles 
that have an axisymmetric metric but a not axisymmet- 


4 L 2 is the Lebesgue space of squared integrable functions. 
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ric embedding. While we can predict the magnitude of 
the principal curvatures, the principal directions of cur¬ 
vature are dictated by imperfections in both experiments 
and numerics when the target metric is hyperbolic. 


C. Variation of Geometric Composition 

We then tested the models (analytical and numerical) 
for other values of R/R e a. When this ratio is between 0 
and 0.5, the Gaussian curvature cannot be approximated 
as homogenous throughout the disk but attains two con¬ 
stant values in the inner disk and the annulus. However, 
the metric does not diverge much from the constant cur¬ 
vature solution, and the analytical model gives a result 
that is in good agreement with the mean Gaussian curva¬ 
ture of the disk, which is exactly what we measured ex¬ 
perimentally. The inset in figure [2] shows the agreement 
between experiments, numerics and analytics as R/R e a 
varies. Notice that the analytical model always overes¬ 
timates the Gaussian curvature since it is based on a 
membrane approximation. The numerical model shows 
that, when the target metric is hyperbolic, wrinkles arise 
below a critical thickness as also shown in[28l We then 
expect our hypothesis of constant Gaussian curvature to 
hold in a finite range of thicknesses. 


III. RESIDUAL SWELLING 
A. Experiments 

Residual swelling is a fairly more complicated phe¬ 
nomenon than the geometrical confinement that dictated 
the shape change in the simplified mechanical problem. 
Swelling-induced deformations cannot be seen as distor¬ 
tions, as they are related to both the elastic properties of 
the gel, and the chemical conditions of the residual free 
polymer chains. Moreover, in this case swelling is driven 
by the concentration gradient of these chains across the 
entire structure. We used the circular molds of the me¬ 
chanical analog to cast samples with polyvinylsiloxane 
as shown in figure [l] (PVS - Zhermack Elite Double 32 
for the annulus and Zhermack Elite Double 8 for the in¬ 
ner disk). Both elastomers are incompressible (y = 0.5) 
and their Young’s moduli were measured as 0.96 MPa 
and 0.23 MPa for PVS 32 (annulus) and PVS 8 (inner 
disk), respectively. The inner disks (radius R) and the 
annuli (radii R and R e ) were geometrically compatible 
so that they could be bonded without pre-stretch. Once 
released from the molds, the geometric composites were 
flat, and the plates morphed into curved disks over time 
due to residual swelling - the flow of free polymer chains 
from high density regions (softer gel, disk) to low density 
regions (stiffer gel, annulus). To study the influence of 
R/R e on the morphing process, we fixed the radius of the 
whole disk to 12 mm and varied the radius of the inner 




r/R e a 

FIG. 2: (a) A schematic of the metric of a circular disk in po¬ 
lar coordinates changing by a as the disk is stretched homoth- 
etically. (b) A schematic of the relevant radii, (c) Schematics 
of the metrics of an annular ring and a circular disk stretched 
independently, along with images of the resulting shell denot¬ 
ing the principal directions of curvature and Gaussian nor¬ 
mal coordinates, (d) A plot of the realized metric (analyti¬ 
cal as solid black and grey curves, numerical as dashed red 
curves, experiments as black symbols) normalized by R 2 a 2 vs. 
the normalized radial coordinate for prescribed target metrics 
(blue curves) when R/R e a — 0.7, which result in positively 
and negatively curved shells. The inset shows the normalized 
Gaussian curvature vs. the radii ratios (analytical as solid 
black and grey lines, numerical as red symbols, experiments 
as black symbols). 


disk from 5 mm to 11 mm casting 7 disks with differ¬ 
ent R/R e . We measured the time evolution of the Gaus¬ 
sian curvature of each disk with the same procedure used 
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for the mechanical analog, repeated every three hours. 


B. Residual Swelling of Geometric Composites 


While this problem couples nonlinear geometric me¬ 
chanics with elastomer swelling, we can provide insight 
into this process by incorporating swelling into our me¬ 
chanical analogy. The stretching ratio a now dictates 
the metric that each part of the disk would realize upon 
swelling if it were free (not bonded to the other). The 
inner disk and the annulus would like to shrink and swell, 
respectively, as molecules are flowing from the former to 
the latter. We assume that if the annulus would like to 
swell by a factor <a, the inner disk would like to shrink 
by a factor a ~ x . Incorporating the difference between 
the two Young’s moduli, the functional in equation © is 
modified as 


/' 


u.*i ) 

a 2 r" E d JJ J 


f 


(aee - a 2 r 2 ) 2 


dr. 


( 5 ) 

Notice that, since no pre-stretch is applied, the radius 
of the disk is R e , i.e. it coincides with the outer radius 
of the green annulus. The Young’s moduli of the green 
annulus and the pink disk are denoted as E a and E re¬ 
spectively . 5 To analytically determine how a should vary 
with R/R e , we denoted as Cd and c a the concentrations 
of the diffusive species in the disk and in the annulus, 
respectively. Since molecules flowed from the disk to the 
annulus, we fixed c a < Cd and imposed the conservation 
of mass as c eq 7rR 2 = CdirR 2 + c a 7r(R 2 — R 2 ), where c eq 
denotes the concentration at equilibrium. Then, we rea¬ 
soned that the stretching ratio a will be proportional to 
the cubic root of the mass uptake inside the annulus so 
that a 3 — l ~ (c eq — c a )7r(R 2 — R 2 ). Finally, by expressing 
c eq from the mass conservation, we got 


a = 


1 + V (Cd - Ca) 





where 77 is a proportionality coefficient having the dimen¬ 
sion of the inverse of a concentration and representing 
the link between mass uptake and stretch . 6 The pres¬ 
ence of a concentration gradient of polymer chains with 
a polydisperse molecular weight makes identifying this 
parameter difficult, and beyond the scope of this work. 
Qualitatively, the bigger the free chains, the higher 77 
should be. Notice that a is equal to 1 when the mass up¬ 
take is zero, that is when the structure is homogeneous 


5 Their ratio is roughly equal to 4 and its variation with swelling 
is neglected. 

6 Similar to the hydrophilicity coefficient introduced in t30l to de¬ 
scribe stretching induced by cation’s motion in ionic polymer- 
metal composites. 


(R/R e = 0 or 1 ). This is the important difference with 
the mechanical analog, since 77 is not known for swelling 
as it depends on the material and chemical properties of 
the elastomers; we therefore used it as a fitting param¬ 
eter. Figure [3] shows the stationary values of the Gaus¬ 
sian curvature of the seven disks after residual swelling 
as obtained in the experiments (triangles), numerics (cir¬ 
cles) and analytics (solid curve). Numerical and analyt¬ 
ical results are obtained by using (| 6 | in © and setting 
r](cd — c a ) equal to 0.54 that sets <a max = 1.043 from ©• 
The three linear regimes identified in figure [3] point out 
that the maximum of the Gaussian curvature is not at¬ 
tained for R/R e = 1/2 but for R/R e ~ 0.77. A similar 
result was obtained for the mechanical analog as can be 
seen in the inset in figure [2] These two observations let 
us conclude that both the dimensionality of the swelling 
and simple geometry shift the maximum of the Gaussian 
curvature to high radii ratios instead of R/R e = 1 / 2 . 
By the conservation of mass, it can be demonstrated 
that if the swelling had been ID, the maximum mass 
uptake would have been attained at R/R e = 1 / 2 ; if it 
had been 3D, the maximum would have been attained 
at R/R e = 1 / 2 1 / 3 ~ 0 . 8 . In our 2 D case, the maximum 
is attained when R/R e = l/y/2 ~ 0.7 as also experi¬ 
ments showed. So, in general, if n is the dimensionality 
of the swelling, the maximum mass uptake is achieved 
for R/R e = 1 / 2 1 / 72 . The agreement among experiments, 
numerics and analytics is quite good, and it is remark¬ 
able that the analytical model captures the linear regimes 
with the same slopes. The closed form analytical solution 
of the problem is cumbersome, but it may be simplified 
by noticing that <a max — 1 , which allows us to perform 
a Taylor expansion in terms of <a max — 1. At the leading 
order, defining E = E a /Ed and R = R/R e , it reads 


KRl 


96(1 


®max)^ 


(l- R 2 )(l- R 3 ) 
R 6 (l — E) +E 


( 7 ) 


which we think can be interpreted as the 2D analog of 
Timoshenko’s formula for beams f25] as it expresses how 
the dimensionless Gaussian curvature varies with mate¬ 
rial and geometric ratios (see Appendix). To the best of 
our knowledge, this is the first analytical formula relat¬ 
ing the Gaussian curvature of a geometric composite to 
its material and geometrical properties, i.e. moduli and 
radii ratios. This simplified expression is represented in 
figure [3] as a grey dashed line, and is very close to the 
full solution: the linear regimes highlighted in the fig¬ 
ure are in excellent agreement with our Timoshenko-like 
formula. The physical interpretation of our first order 
Taylor approximation is that the strains are assumed to 
be small, which is the same limit that Timoshenko ob¬ 
tained his formula within. It is worth noting that, unlike 
thermal stretches in uniformly heated bimetallic strips, 
the stretching ratio a should depend on the elastic prop¬ 
erties of the geometric composites, as discussed in I3T1 
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r/R e r/R e r/R e r/R e 

FIG. 3: (a) Stationary Gaussian curvature vs R/R e achieved 
by residual swelling. Three different regimes may be identi¬ 
fied as the radii ratio varies. In the first one (I) the curvature 
is very small. The second one (II, linear increasing ) shows a 
linear scaling of the Gaussian curvature with the radii ratio 
up to its maximum when the third regime (III, linear de¬ 
creasing) starts with a steeper linear scaling of the curvature 
now decreasing to zero. The solid black line is the analyti¬ 
cal solution, the dashed grey line is its Taylor approximation 
(eq. 0 ); triangles and circles represent experimental and nu¬ 
merical results, respectively, (b) Deformed shapes for four 
different radii ratios and corresponding experimental profiles 
of the Gaussian curvature. 


C. Swelling Dynamics 

Our model successfully captures the steady-state mor¬ 
phology of residually swollen plates. Unlike the mechani¬ 
cal analog presented earlier, the residual swelling process 
adds a time-dependency to the deformation. The exper¬ 
imental results in figure [4] (a) show the time evolution 
of the Gaussian curvature of disks with seven different 
R/R e , and the shape evolution contains two notable fea¬ 
tures: 1.) there is a critical activation time , i.e. the time 
it takes for the structure to start deforming that depends 
on R/R e , and 2.) following actuation, the disks deform 
in a diffusive manner. 

We assume that the swelling dynamics may be de¬ 
scribed as a diffusive process with a Fourier-like differen¬ 


tial equation. The main features of a diffusive equation 
are that it is of the first order in time, giving rise to 
transients that are described by the exponential of time 
up to the steady-state, and of the second order in space 
(quasi-ID in our case). As we are studying the transient 
by looking at a homogeneous field - the Gaussian curva¬ 
ture K - we focus on its variation with time, and note 
that the dashed lines in figure [4] (a) correspond to an 
exponential of time (K ste ^ y R 2 (l — e _t / r )), as expected. 
Figure [4] (b) shows that the activation time varies with 
R/R e as t a ~ A e -BR/R e ^ w h ere a an q g are positive 
real numbers equal to 21038 h and 10.862, respectively, 
in our case. 7 This numerical fitting is shown in the plot 
as a straight dashed grey line. Following activation, the 
diffusive shape change is characterized by the time con¬ 
stant r that dictates the time scale of the transient as 
it is the time at which the Gaussian curvature reaches 
the 63% of its stationary value. Figure [4] (a) shows that 
the disk corresponding to R/R e = 11/12 is faster than 
the other geometries; indeed, we found that its time con¬ 
stant is approximately r ^ 40 h whereas all the other 
disks have r ^ 90 h as shown in figure [4] (c). 

As this is a diffusive process, we expect the dynam¬ 
ics will scale with the square of the characteristic length 
scale in the problem. We believe the observed difference 
in dynamics is the result of a change in the relevant char¬ 
acteristic length scale, i.e. from the total radius of the 
disk to the width of the annulus. The characteristic time 
scale (time constant) of a diffusive process is equal to 
r = £ 2 /D, where £ is the characteristic length and D is 
the diffusivity. While the latter is a property of the mate¬ 
rials, the former is dictated by geometry. To identify the 
characteristic length, we examined a simpler geometry - 
a bilayer beam. Using the same materials as in section 
3.1, we prepared a bilayered beam of equivalent thick¬ 
ness (figure [4] a - inset), and measured the time it took 
equilibrate into an arch, finding r^ eam ^ 10 h represented 
by a diamond in figure [4] (c). In this case, a reasonable 
assumption for the characteristic length is £ ~ h, where 
h is the total thickness of the beam as shown in the inset 
of figure [4] (a). Since the materials for the beam and disk 
are the same, they share the same value of D. There¬ 
fore, we compared the ID diffusion in the beam with the 
2D diffusion in the disk with R/R e = 1/2 and assumed 
£ ~ R e obtaining Tdi s k = (i? e /^) 2 7beam — 90 h: this ana¬ 
lytical estimate is shown in figure [4] (c) as a square and 
it excellently predicts the experimental time constant of 
the disk. The experimental data show a decay of the 
time constant as the radii ratio approaches 1 that we in¬ 
terpret as the result of a decaying characteristic length, 
which represents the portion of the radius where swelling 
is actually taking place. When R/R e ~ 0.5, our approx¬ 
imation £ rsj R e is a good estimate for the characteristic 


7 Error bars correspond to ±3 h since we measured K every three 
hours. 
























length but when R/R e —>> 1 the inner area of the disk 
is shielded from swelling and the characteristic length is 
smaller than R e . We suggest that, as R/R e 1, the 
characteristic length approaches the width of the annu¬ 
lus. 


a. 




R/Re R/Re 


FIG. 4: (a) Normalized Gaussian curvature versus time for 
the different radii ratios, (b) Activation times versus R/R e 
(triangles) and numerical fitting (dashed curve), (c) Time 
constants of the disks versus R/R e (circles), time constant 
of the beam (diamond) and analytical estimate of the time 
constant of the disk corresponding to R/R e = 1/2 (square) 
assuming £ ~ R e . 


IV. CONCLUSIONS 


how the Gaussian curvature is dictated by geometry. The 
strength of the model is indeed its analytical tractabil- 
ity that results from the assumption of a homogeneous 
Gaussian curvature throughout the disk. The agreement 
among experiments, numerics, and analytics is excellent 
even when the Gaussian curvature is not homogeneous 
because the analytical model provides a mean Gaussian 
curvature as a result, which is important for the design 
of actuators. 

We then employed the analytical model of the mechan¬ 
ical analog to study the morphing of geometric compos¬ 
ites by approximating the swelling as a distortion. By us¬ 
ing the conservation of mass, we analytically determined 
how the mass uptake should vary with the radii ratio, 
and assumed a linear proportionality between the mass 
uptake and the volume variation. The agreement among 
experiments, numerics and analytics is quite good and 
each approach identified three regimes for the Gaussian 
curvature as a function of the radii ratio: it is remark¬ 
able that the model catches these regimes and their linear 
features. Then, by assuming small stretches (a ~ 1), we 
simplified the cumbersome analytical solution and pro¬ 
vided the first 2D extension of the Timoshenko’s formula 
for beams. Finally, we studied the swelling dynamics and 
identified two different characteristic lengths depending 
on geometry. 

We think that the proposed model improved the un¬ 
derstanding of the complex interplay among geometry, 
mechanics, and swelling. Additionally, the experiments 
demonstrate a robust and scalable means for the growth- 
actuated manufacturing of elastic shells - a material that 
is traditionally difficult to prepare via additive and reduc¬ 
tive manufacturing techniques. It is important to note 
that while residual fluid within the crosslinked elastomer 
drives the diffusion and swelling of the structure, the ma¬ 
terial behaves like an elastic solid, rather than a swollen 
gel. We expect this experimental procedure to translate 
to any combination of material-compatible elastomers 
where a gradient of small molecules can be programmat¬ 
ically prescribed. This may provide the foundation for 
an inkjet-like approach to 3D printing whereby small 
molecule fluids can be locally applied to a flat elastic 
sheet, allowing controlled diffusion to dictate the result¬ 
ing growth pattern. Careful selection of the initial geom¬ 
etry will allow this technique to be used for generating 
regions of high curvature - or folds - which may form the 
basic building blocks for the growth of origami structures. 


We have studied the morphing of geometric composites 
from flat plates into curved, three-dimensional shapes. 
The geometric composites morph by residual swelling, 
a phenomenon that we gained insight into by consider¬ 
ing a mechanical analog that copies their geometry and 
morphs by geometrical confinement. The morphing prob¬ 
lem of the mechanical analog is purely geometrical, and 
we developed an analytical model following the theory 
of non-Euclidean plates, which quantitatively describes 
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Appendix A: Derivation of the stretching energy of 
the mechanical analog 

Equation 0 is derived from 0 by computing the 
traces of the tensors a — a and (a — a) 2 using their po¬ 
lar coordinates. By definition, the trace of a tensor is 
the result of the contraction between the tensor and the 
metric tensor. As an example, the trace of a — a may be 
computed in polar coordinates as 

tr (a - a) = a a0 (a a/3 - a a p), 

where d a/3 denotes the contravariant components of the 
target metric tensor that may be thought as the entries 
of the inverse of the matrix [a a p]- So, a rr = / _2 (r), 
a r0 = 0 and a 60 = f~ 2 {r)r~ 2 . Simple algebra yields: 

tr (a-a) = a aP a a0 - 2 = f~ 2 (r)a rr + f~ 2 (r)r~ 2 a e o - 2 . 

Then, if CLrprp (l rr ^ we get 

tr 2 (a — a) = ( f~ 2 (r)r~ 2 a ge - l) 2 = (^ - l) . 

\aee ) 

In our case, the trace of (a —a) 2 turns out to be equal to 
tr 2 (a — a) so that 

tr 2 (a — a) + tr (a — a) 2 = 2 . 

\aoo ) 

Finally, to compute the stretching energy, we have to 
evaluate that is 

VW\ = f 2 {r)r- 

By using the last two formulas, equation 0 is recovered 
after some rearrangements, discarding a factor 2 since we 
are just interested in the minimization of the energy. 


Notice that, since we are using Gaussian normal coordi¬ 
nates, a pp = a rr (dr/dp) 2 a rr f~ 2 (r ) = 1 and clqq is 
the only metric coefficient contributing to the Gaussian 
curvature. 

Appendix B: Timoshenko—like formula 

Timoshenko wrote his formula for the bending of bi¬ 
metal thermostats as[25] 

^ 6A£(l + ra) 2 

3(1 + m) 2 + (1 + mn)(m 2 + ’ 

where k is the curvature of the beam, h is its thickness, 
Ae is the difference between the thermal strain in the two 
layers, m is the ratio between the thicknesses of the two 
layers and n is the Young moduli’s ratio. Equation 0 
resembles Timoshenko’s formula in that it expresses a di¬ 
mensionless curvature, albeit a Gaussian curvature, as a 
function of geometric and material ratios. Notice that the 
characteristic length is the one along which the compos¬ 
ite is non homogeneous: the total thickness in the beam 
and the total radius in the disk, which may be thought 
of being obtained by rolling the bilayer beam. 


1. Gaussian normal coordinates 


As a rr ~ a rr , the Gaussian normal coordinate p is 
computed as 



r, r < R 

R + a(r — R ), r > R 


2. Gauss’s Theorem 


Since the polar coordinates are orthogonal, the Gauss’s 
Theorem may be written as 


K = - 


1 

y/CL pp dQQ 



,p 
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